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Abstract. We apply the semidefinite programming method to derive 
bounds for projective codes over a finite field. 

1. Introduction 

In network coding theory, as introduced in [1] , information is transmitted 
through a directed graph. In general this graph has several sources, sev- 
eral receivers, and a certain number of intermediate nodes. Information is 
modeled as vectors of fixed length over a finite field ¥ q , called packets. To 
improve the performance of the communication, intermediate nodes should 
forward random linear F^-combinations of the packets they receive. This is 
the approach taken in the non-coherent communication case, that is, when 
the structure of the network is not known a priori [13]. Hence, the vector 
space spanned by the packets injected at the source is globally preserved in 
the network when no error occurs. This observation led Koetter and Kschis- 
chang [15] to model network codes as subsets of projective space V(F™), the 
set of linear subspaces of F™, or of Grassmann space G q (n,k), the subset 
of those subspaces of IF™ having dimension k. Subsets of V(¥^) are called 
projective codes while subsets of the Grassmann space will be referred to as 
constant- dimension codes or Grassmann codes. 

As usual in coding theory, in order to protect the system from errors, it 
is desirable to select the elements of the code so that they are pairwise as 
far as possible with respect to a suitable distance. The subspace distance 
between U and V 

d s (U, V) = dim(C7 + V) - dim(C7 n V) = dim U + dim V - 2 dim(U D V) 
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was introduced in [15J for this purpose. It is natural to ask how large a code 
with a given minimal distance can be. Formally, we define 

A q (n,d) := max{|C[ : C C P(FJ), d s (C) > d} 
A q (n, k,2S) := max{|C| : CcG q (n,k), d s {C)>25} 

where ds{C) denotes the minimal subspace distance among distinct elements 
of a code C. In this paper we will discuss and prove upper bounds for A q (n, d) 
and A q (n, k, 28). 

1.1. Bounds for A q (n, k, 25). Grassmann space Q q (n,k) is a homogeneous 
space under the action of the linear group GL n (¥ q ). Moreover, the group 
acts distance transitively when we use the subspace distance; the orbits of 
GL n (¥ q ) acting on pairs (U, V) of Gq(n, k) are characterized by the subspace 
distance d(U,V). In other words, Grassmann space is two-point homoge- 
neous under this action. 

Due to this property, codes and designs in Q q (n, k) can be analyzed in the 
framework of Delsarte's theory, in the same way as other classical spaces in 
coding theory, such as Hamming space and binary Johnson space. In fact, 
Gq(n, k) is a q-analog of binary Johnson space; see [7J. The linear group plays 
the role of the symmetric group for the Johnson space, while the dimension 
replaces the weight function. 

The classical bounds (anticode, Hamming, Johnson, Singleton) have been 
derived for the Grassmann codes [15, 26, [27]. The more sophisticated Del- 
sarte linear programming bound was obtained in [7j. However, numerical 
computations indicate that it is not better than the anticode bound. More- 
over, the Singleton and anticode bounds have the same asymptotic behavior 
which is attained by a family of Reed-Solomon-like codes constructed in |15j 
and closely related to the rank-metric Gabidulin codes. 

1.2. Bounds for A q (n, d). In contrast to G q (n, k), the projective space has a 
much nastier behavior, essentially because it is not two-point homogeneous. 
In fact it is not even homogeneous under the action of a group. For example, 
the size of balls in this space depends not only on their radius, but also on 
the dimension of their center. Consequently, bounds for projective codes are 
much harder to obtain. Etzion and Vardy in [10] provide a bound in the 
form of the optimal value of a linear program, which is derived by elementary 
reasoning involving packing issues. Up to now the Etzion- Vardy bound is 
the only successful generalization of the classical bounds to projective space. 

In this paper we derive semidefinite programming bounds for projective 
codes and compare them with the above mentioned bounds. 

In convex optimization, semidefinite programs generalize linear programs 
and one can solve them by efficient algorithms [23], [24]. They have nu- 
merous applications in combinatorial optimization. The earliest is due to 
Lovasz [18] who found a semidefinite programming upper bound, the theta 
number, for the independence number of a graph. 
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Because a code with given minimal distance can be viewed as an indepen- 
dent set in a certain graph, the theta number also applies to coding theory. 
However, because the underlying graph is built on the space under consider- 
ation, its size grows exponentially with the parameters of the codes. So by 
itself the theta number is not an appropriate tool, unless the symmetries of 
the space are taken into account. A general framework for symmetry reduc- 
tion techniques of semidefinite programs is provided in [3]. For the classical 
spaces of coding theory, after symmetry reduction, the theta number turns 
out to be essentially equal to the celebrated Delsarte linear programming 
bound. For projective spaces, the symmetry reduction was announced in 
[5] (see also [2]). The program remains a semidefinite program (it does not 
collapse to a linear program) but fortunately it has polynomial size in the 
dimension n. 

The relationship between Delsarte's linear programming bound and the 
theta number was recognized long ago in [21] and [20]. Recently, more appli- 
cations of semidefinite programming to coding theory have been developed, 
see [22], [I], [25], [H] and the survey [2]. 

1.3. Organization of the paper. In Section [2] we review the classical 
bounds for Grassmann codes and the Etzion-Vardy bound for projective 
codes. In Section [3] we present the semidefinite programming method in 
connection with the theta number. We show that most of the bounds for 
Grassmann codes can be derived from this method. In Section [J] we reduce 
the semidefinite program by the action of the group GL n (¥ q ). In Section [5] 
we present numerical results obtained with this method and we compare 
them with the Etzion-Vardy method for q = 2 and n < 16. Another distance 
of interest on projective space, the injection distance, was introduced in [T7J. 
We show how to modify the Etzion-Vardy bound as well as the semidefinite 
programming bound for this. 

2. Elementary bounds for Grassmann and projective codes 

2.1. Bounds for Grassmann codes. In this section we review the classi- 
cal bounds for A q (n, k,25). We note that the subspace distance takes only 
even values on the Grassmann space and that one can restrict to k < n/2 
by the relation A q (n, k, 25) = A q (n, n — k, 25), which follows by considering 
orthogonal subspaces. 

We recall the definition of the g-analog of the binomial coefficient that 
counts the number of ^-dimensional subspaces of a fixed n-dimensional space 
over ¥ q , i.e. the number of elements of Q q {n, k). 

Definition 2.1. The g-ary binomial coefficient is defined by 



n 




{q n - 1) . . . ( g "- fc + 1 - 
(q k - 1) . . . (q - 1) 
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2.1.1. Sphere-packing bound. 

Q q (n,k)\ 



(1) 



A q (n,k,25) < 



La 



\ Bk (5-i)\ Yl% im [L] q r^r 2 

It follows from the well-known observation that balls of radius 5 — 1 centered 
at elements of a code C C Q q (n,k) with minimal distance 25 are pairwise 

disjoint and have the same cardinality Yll^=o^ 2 ^ [m] ["ml ^ ■ 



2.1.2. Singleton bound (15] . 

(2) A q (n,k,25)< 



n-S + 1 
k-5 + l 



It is obtained by the introduction of a "puncturing" operation on the code. 

2.1.3. Anticode bound [26\ . An anticode of diameter e is a subset of a metric 
space whose pairwise distinct elements are at distance less or equal than e. 
The general anticode bound (see [6]) states that, given a metric space X 
which is homogeneous under the action of a group G, for every code C C X 
with minimal distance d and for every anticode A of diameter d— 1, we have 



\C\< 



\x\ 



Spheres of given radius r are anticodes of diameter 2r. So if we take A to 
be a sphere of radius 5— 1 in Q q (n, k) we recover the sphere-packing bound. 
Obviously, to obtain the strongest bound, we have to choose the largest 
anticodes of given diameter, which in our case are not spheres. Indeed, the 
set of all elements of Q q (n, k) which contain a fixed (k — 5 + l)-dimensional 
subspace is an anticode of diameter 25 — 2 with [""^f" 1 ]^ elements and 
which is in general larger than the sphere of radius 5 — 1. Moreover Frankl 
and Wilson proved in [TI] that these anticodes have the largest possible size. 
Taking such A in the general anticode bound, we recover the (best) anticode 
bound for G q (n, k): 



A q (n,k,26) < 



["■] 



(3) 



n— k+S— 1] 
5-1 J 



(q n -l)(q 



n-1 



!)...(? 



n— k+S 



(q k ~ 1)(« 



fc-i 



l)...(q*-l) 



It follows from the previous discussion that the anticode bound improves 
the sphere-packing bound. Moreover, the anticode bound is usually stronger 
than the Singleton bound, with equality only in the cases n = k or 5 = 1, 
see [27] . 
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2.1.4. First and second Johnson-type bound 

(a n - l)(a k 

(4) A q (n,k,2S)< 1 Mi 
as long as (q k — l) 2 

(5) A q (n,k,25)< 



(q k - l) 2 - (q n - l)(q k - 5 - 1) 

(q n - l){q k - & - 1) > 0, and 
n _ i 



1 



AJn-l,k-l,25) 



These bounds were obtained in [27] through the construction of a binary 
constant-weight code associated to every constant-dimension code. Iterating 
the latter, one obtains 



(6) A q (n,k,28)< 



1 



_,n— 1 



1 



1 



7 n— k+8 



1 



If the floors are removed from the right hand side of Q, the anticode bound 
is recovered, so © is stronger. In the particular case of S = k and if 
n ^ mod k, ([6]) was sharpened in [10] to 

- 1 



(7) 



A q (n,k,2k) < 



1. 



For 5 = k and if k divides n, we have equality in ([6]), because of the existence 
of spreads (see [10] ) 

A q (n,k,2k) 



1 



<f - 1 

Summing up, the strongest upper bound for constant dimension codes re- 
viewed so far comes by putting together ([6]) and ([7]): 



Theorem 2.2. If n — k ^ mod 5, then 

7 n— fc+<5+l i 

A q {n,k,25)< " 
otherwise 



q n - 1 




_q k - 1 





7 5+l _ i 



-.n—k+8 



A q (n,k,28) < 





- 1 




qk 


- 1 





n-k+8+l 



7 5+l 



1 



«— k+8 



1 



2.2. A bound for projective codes. Here we turn our attention to pro- 
jective codes whose codewords have not necessarily the same dimension, and 
we review the bound obtained by Etzion and Vardy in [10]. The idea is to 
split a code C into subcodes Ck = C n G q (n,k) of constant dimension, and 
then to derive linear constraints on the cardinality |C&|, coming from packing 
constraints. 

Let B(V, e) := {U G P(F") : d s (U, V) < e} denote the ball with center V 
and radius e. If dim!/ = i, we have 



\B(V,e) 



EE 

1=0 j=0 



i 




n — i 


j_ 


q 


A- 3. 
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We define c(i, k,e) := \B(V,e) H Q q (n,k)\ for V of dimension i. It is not 
difficult to prove that 

min{/c,i} 

(8) c(i,k,e)= Yl 



i 




n — i 


j. 




k - j 



q(i-j)(k-j) _ 

1 



Theorem 2.3 (Linear programming bound for codes in 'P(Fg), |10j). 

n 

A g (n, 2e + 1) < max | Xk : Xk < fc, 2e + 2) /or fc = 0, . . . , n, 

fc=o 

n r- -i 

c(i, fc, e)xj < /or fc = 0, . . . , n j 
i=o L J 9 

Proof. For C C P(F~) of minimal distance 2e + 1, and for fc = 0, . . . , n, we 
introduce = \C H £7 9 (n, fc)|. Then X}fc=o x k = 1^1 anc ^ eacn x fc represents 
the cardinality of a subcode of C of constant dimension fc, so it is upper 
bounded by A q {n,k,2e + 2). Moreover, balls of radius e centered at the 
codewords are pairwise disjoint, so the sets B(V, e) n G q (n, fc) for V £ C are 
pairwise disjoint subsets of G q (n, fc). So 

Y \B(V,e)nG q (n,k)\ < \G q (n,k)\. 
v&c 

Because \B(V, e) D fc)| = c(i, fc, e) if dim(V) = i we obtain the second 
constraint 

n i- -i 

y^c(i,fc,e)xi < . 

i=0 L J g 

So |C| is at most the optimal value of the linear program above. □ 

Remark 2.4. Of course, in view of explicit computations, if the exact value 
of A q (n,k,2e + 2) is not available, it can be replaced in the linear program 
of Theorem \2.3\ by an upper bound. 

3. The semidefinite programming method 

3.1. Semidefinite programs. A (real) semidefinite program is an opti- 
mization problem of the form: 

sup{ (C,Y) : YhO, (Ai,Y) = bi for i = I,..., m }, 

where 

C, Ai, . . . , A m are given real symmetric matrices, 
&1, . . . , b m are given real values, 

Y is a real symmetric matrix, which is the optimization variable, 
(A, B) = trace(^4-B) is the inner product between symmetric matri- 
ces, 

Y y denotes that Y is symmetric and positive semidefinite. 
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This formulation includes linear programming as a special case when all 
matrices involved are diagonal matrices. When the input data satisfies some 
technical assumptions (which are fulfilled for our application) then there are 
polynomial time algorithms which determine an approximate optimal value. 
We refer to [23] and [24] for further details. 

3.2. Lovasz' theta number. In [18], Lovasz gave an upper bound on the 
independence number a(Q) of a graph Q = (V, E) as the optimal value $(0) 
of a semidefinite program: 

Theorem 3.1 ([18]). 

a(Q) < #(Q) ■= max{ £ F(x,y) : F G R VxV , F y 0, 
(9) J2F(x,x) = l, 

F(x,y) = 0ifxy£E} 

Here we can write max instead of sup because one can show using duality 
theory of semidefinite programming that the supremum is attained: the 
Slater condition [24, Theorem 3.1] is fulfilled. 

In the above and all along this paper, we identify a matrix indexed by a 
given finite set V with a function defined on V 2 . The program given in ([2]) 
is one of the many equivalent formulations of Lovasz' original "&(Q)- If the 
constraint that F only attains nonnegative values is added, the optimal value 
gives a sharper bound for a(Q). Traditionally this semidefinite program is 
denoted by |2T]. 

This method applies to bound the maximal cardinality A(X, d) of codes 
in a metric space X with prescribed minimal distance d. Indeed A(X, d) = 
a(Q) where Q is the graph with vertex set X and edges set {xy : < 
d(x,y) < d}. So, we obtain: 

Corollary 3.2 (The semidefinite programming bound). 
(10) 

A{X,d)<max\ £ F(x,y) : F G R XxX , F y 0, F > 0, 

E F{x,x) = 1, 
xex 

F(x,y) = if0<d(x,y) <d} 

= min|i/A : F £ R XxX , F - X y 0, 
(11) ^ F(x,x)<t forallxeX, 

F{x,y)<0 ifd{x,y)>d} 

The second semidefinite program (jllj) is the dual of (jlOp . Furthermore, 
by weak duality, any feasible solution of the semidefinite program in (jllj) 
leads to an upper bound for A(X, d). 
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3.3. Bounds for Grassmann codes. In two-point homogeneous spaces, 
the semidefinite program in (llli collapses to the linear program of Delsarte, 
first introduced in [6] in the framework of association schemes. This fact 
was first recognized in the case of the Hamming space, independently in [20 
and [21] . We refer to [3] for a general discussion on how (fTOj) and (fTTj) reduce 
to linear programs in the case of two-point homogeneous spaces. 

Grassmann space G q (n, k) is two-point homogeneous for the action of the 
group G = GL n (¥ q ) and its associated zonal polynomials are computed in 
[7]. They belong to the family of g-Hahn polynomials, which are g-analogs 
of the Hahn polynomials related to the binary Johnson space. 

Definition 3.3. The q-Hahn polynomials associated to the parameters n, s, t 
with < s < t < n are the polynomials Qe(n, s, t; u) with < i < min(s, n — 
t) uniquely determined by the properties: 

(a) Qi has degree t in the variable [u] = q^ u ["] 

(b) They are orthogonal polynomials for the weights 

< i < min(s, n — t) w(n, s, t; i) = 

(c) Q e (0) = 1. 

To be more precise, in the Grassmann space Q q (n,k), the zonal polyno- 
mials are associated to the parameters s = t = k. The other parameters will 
come into play when we analyze the full projective space in Section HI The 
resulting linear programming bound is explicitly stated in [7j: 

Theorem 3.4 (Delsarte's linear programming bound [7J). 

A q (n,k,25) <min{l + /i + ■■■ + /* : /« > for i = 1, . . . , k, 

F(u) <0 foru = 6,...,k\, 

where F(u) = 1 + Ym=i fiQi{ u ) an d Qi( u ) = Qi(n, k, k;u) as in Defini- 
tion\3M 

In order to show the power of the semidefinite programming bound, we 
will verify that most of the bounds in Section [2] for Grassmann codes can 
be obtained from Corollary 13.21 or Theorem 13.41 In each case we construct 
a suitable feasible solution of (fTTI) . 

3.3.1. Singleton bound. We fix an arbitrary subspace w of F™ of dimension 
n — 5 + 1. We consider a function </> : Q q (n,k) —> {u (Z w : dim(n) = k — 6+1} 
such that <j>(x) C x for all x. Clearly dim(x fl w) > k — S + 1. In the case 
of equality, we set <p{x) = x fl w. If dim(x n w) > k — 8 + 1, (p(x) is chosen 
arbitrarily among the (k — 5 + l)-dimensional subspaces of x fl w. 



n — s 
t — s + i 



i(t—s+i) 
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We define the function 

F(x,y) 



1 if fa) = cb(y) 
otherwise 



£ l(<P(x) = u)l(4>(y) = u) 



uGw 
dim(u)=fc— 5+1 



where \(4>(x) = u) denotes the characteristic function of the set {x £ 
G g (n, k) : 4>(x) = u}. Then, F is obviously positive semidefinite, and (F, i, A) 
is a feasible solution of (1111) where t = 1 and 



A 



-2 

E F ( x >v) 

1 (x,y)eg q (n,k) 2 

E ( E 

dim(M)=fc— 5+1 



l(0(x) =«))' 

xeG q (n,k) 



It follows from Cauchy-Schwarz inequality that A > so the Single- 

ton bound ([21) is recovered from (jlip . 

3.3.2. Sphere-packing and anticode bounds. The sphere-packing bound and 
the anticode bound in Q q (n, k) can also be obtained directly, with 

F{x,y)= ^2 l B (z,5-l){x)lB{z,5-l){y)i 
dim(z)=fc 

and 

F(x,y)= HzCx)l(zCy). 

dim(z)=fc— 5+1 

In general, the anticode bound \C\ < |X|/|.A| can be derived from (ITT]) , using 
the function F(x,y) = E 9 gg 1 A(gx)lA(9y)- 

3.3.3. First Johnson-type bound. We want to apply Delsarte's linear pro- 
gramming bound of Theorem 13.41 with a function F of degree 1, i.e. F(u) = 
foQo(u) + fiQ\{u). According to [7] the first g-Hahn polynomials are 

/ ( Q n - lVl - n~ u ) 
Q (u) = l , Ql (u) ={l- * _ ^ qn J_ ^ 

In order to construct a feasible solution of the linear program, we need 
fo, fi > for which F(u) = fo + fiQi(u) is non-positive for u = 5, . . . , k. 
Then 1 + /1//0 will be an upper bound for A q (n, k, 25). As Qi(u) is decreas- 
ing, the optimal choice of (fo, fi) satisfies F(5) = 0. So fi/fo = —1/Q\(S) 
and we need Qi(5) < 0. We obtain @: 



fa <?l(i) (q k -l) 2 -(q"-l)(q>- l -iy 
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3.3.4. Second Johnson-type bound. Here we find an inequality for the op- 
timal value B q {n,k,25) of the semidefinite program (|lip in the case X = 
Q q {n,k) (with the subspace distance) which resembles ([5]): 

a n - 1 

B q (n, k, 25) < 'L—Bgin - 1, k - 1, 25). 

Let (F, t, A) be an optimal solution for the program (jlip in Q q {n — 1, k — 1) 
relative to the minimal distance 25, i.e. i* 1 satisfies the conditions: -F >z A, 
x) < t, F(x, y) < if d(x, y) > 25, and t/A = - 1, k - 1, 25). We 
consider the function G on Q q {n, k) x £/ g (n, k) given by 

G(x,y)= £ 1(D Cx)l(D Cy)F(xnH D ,yDH D ), 

dim(D)=l 

where, for every one-dimensional space D, Hp is an arbitrary hyperplane 
such that D © Ho = F" It can be verified that the triple (G,t',X r ) is a 
feasible solution of the program (jlip in Q q {n, k) for the minimal distance 25, 

where t' = t[ l [] q and A' = ^[x] g /[i] ? ; thus leading to the upper bound 

+' + a n — 1 a n — 1 
25) < - = = ^B,^- l,k- 1,25). 

Remark 3.5. In [10] another Johnson-type bound is given: 

a n - 1 

^(n, k, 25) < l_ k _ i A q (n - 1, fc, 25), 

which follows easily from the second Johnson-type bound combined with the 
equality A q (n,k,25) = A q (n,n — k,25). Similarly to above, an analogous 
inequality holds for the semidefinite programming bound B q (n,k,25). 

4. Semidefinite programming bounds for projective codes 

In this section we perform a symmetry reduction of the semidefinite pro- 
grams pO]) and (fTT|) in the case of projective space, under the action of the 
group G = GL n (¥ q ). We follow the general method described in [3]. The 
key point is that these semidefinite programs are left invariant under the 
action of G so the set of feasible solutions can be restricted to G-invariant 
functions F. The main work is to compute an explicit description of the 
G-invariant positive semidefinite functions on the projective space. 

4.1. G-invariant positive semidefinite functions on projective spaces 

In order to compute these functions, we use the decomposition of the space 
of real-valued functions under the action of G. We take the following nota- 
tions: 

X = V{¥ n q ), X k = Q q (n,k), R x = {f :X ^R}. 
The space M. x is endowed with the inner product (, ) defined by: 

(f,g) = t^E/^)^)- 
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For k = 0, . . . , n, an element of R Xk := {/ : X k -> R} is identified with the 
element of R^ that takes the same value on X k and the value outside of 
Xk- In this way, we see the spaces R Xk as pairwise orthogonal subspaces of 
R x . 

Delsarte |7) showed that the irreducible decomposition of the R^ fc under 
the action of G is given by the harmonic subspaces H k ^. 

(12) R Xk = H 0:k e H 1)k e ■ ■ ■ e H min{ktn _ k}tk 

Here, H k k is the kernel of the differentiation operator 

5 k : R Xk — ► R Xfc -! 

/ — > [ x -> ^{/(y) : dim(y) =I:,iCi/}] 

and H kt i is the image of H k ,k under the valuation operator 

ip ki : R Xk — ► R x * 

/ — > [ x -> : dim(y) = fe, y C x} ] 

for k < i < n — k. Because is surjective, we have h k := dim(H k k ) = 
[fc]<j ~~ [fe— J ■ Moreover, ip^ commutes with the action of G, so H k i is 
isomorphic to Hk k- Putting together the spaces one gets the global 
picture: 





= R x « C 




2b 


= #o,o £ 


D #o,i 


2i 




#i,i 


x 2 







LfJ 



H. 



o-LfJ 
2 > LfJ 



#0,(71-1) 
#1,(71-1) 



#0/. 



x LfJ - %J.LfJ ® "•• 

Here, the columns give the irreducible decomposition (112p of the spaces R^* . 
The irreducible components which lie in the same row are all isomorphic, 
and together they form the isotypic components 

jn— 2m+l 



T 

-L"n 



Hn 



H 



m,m+l 



H 



m,n—m — 11 m,m 



Starting from this decomposition, one builds the zonal matrices E k (x,y) [H 
Section 3.3] in the following way. We take an isotypic component X k and we 
fix an orthonormal basis (e kkl , ... , e kkhk ) of H k)k . Let e ksi := 4> ks (e kki ). It 
follows from Theorem 3] that (e ks i, . . . ,e ks h k ) is an orthogonal basis of 
H k s and that 



(13) 



(e k sii &ks 



n - 
s 



2k 

- k 



k(s-k) 
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Then we define E k (x,y) G R(»-2**i)x(r»-2*+i) entrywise by 

hie 

(14) E kst (x, y) eksi{x)e kti (y). 

i=l 

We note that [31 Theorem 3.3] requires orthonormal basis in every subspace, 
while in the definition (|14p of E kst we do not normalize the vectors e ks i. 
Because the norms (fl~3|) do not depend on i, but only on k, s, the matrix 
(E' k (x, y)) s ,t associated to the normalized basis is obtained from (E k (x, y)) s ,t 
by left and right multiplication by a diagonal matrix. So the characterization 
of the G-invariant positive semidefinite functions given in [3j Theorem 3.3] 
holds aswell with (|14p : 

Theorem 4.1. F G M XxX is positive semidefinite and G -invariant if and 
only if it can be written as 

\n/2\ 

(15) F(x,y)= Y,( F k,E k (x,y)) 

k=0 

where F k G ]j( n - 2fc +l)x(n,-2fc+l) and Fq, . . . , F\ n / 2 \ are positive semidefinite. 

Now we compute the E k 7 s explicitly. They are zonal matrices: in other 
words, for all k < s, t < n - k, for all g G G,E kst (x,y) = E kst (gx,gy). 
This means that E ks t is a function of the variables which parametrize the 
orbits of G on X x X. It is easy to see that the orbit of the pair (x, y) is 
characterized by the triple (dim(x), dim(y), dim(x n y)). 

The next theorem gives an explicit expression of E ks t, in terms of the 
polynomials Q k of Definition 13.31 

Theorem 4.2. Ifk<s<t<n — k, dim(x) = s, dim(y) = t, 

E kst (x, y) = \X\h k S ~ k ' g *~* q q k{t ' k) Qk{n, s,t;s- dim(x n y)) 

Vtlqlslq 

//dim(x) ^ s or dim(y) ^ t, E kst (x,y) = 0. 

We note that the weights involved in the orthogonality relations of the 
polynomials Q k have a combinatorial meaning: 

Lemma 4.3 ([8]). Given x G X s , the number of elements y G Xt such that 
dim(x Hy) = s — i is equal to w(n, s, t; i). 

Proof of Theorem \4-%\ By construction, E kst (x,y) ^ only if dim(x) = s 
and dim(y) = t, so in this case E kst is a function of (s — dim(x n y)). 
Accordingly, for k < s < t < n — k, we introduce P ks t such that E kst (x, y) = 
P ks t{s — dim(x Dy)) . Now we want to relate P ks t to the g-Hahn polynomials. 
We start with two lemmas: one obtains the orthogonality relations satisfied 
by P ks t and the other computes P ks t(0). 
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Lemma 4.4. With the above notations, 

rt— fe] rn— 2k] 

(i6) Pkst (o) = ^r^v *-**. 



Proof. We have Pkstifi) = E kst {x, y) for all x, y with dim(x) = s, dim(y) = t, 
x C y. Hence, 



Pkst(0) = r«1 -TK Y E kst{x,y) 

Itlqlslq dim(x)=s 



dim(y)=t 
xCy 



-JFT Y ^2 e ksi(x)e kti (y) 



lt\q LsJ q dim(x)=s i=l 
dim(y)=t 
xCy 



h k 

\"] H ^ 

It 1 q LsJ q i=l 

h k 



1 / \ 

PksM = rra1 ni Y Y { Y e ksi{x)jekti(y) 



q LsJ q i=l dim(y)=t dim(x)=s 
xCy 



~m~Y Y ^s,t{ e ksi){y)e k ti{y)- 

LtJ q LsJ q i=l dim(y)=t 



With the relation il) st o = [g_fc] V'fct) 



^st(efesi) = V'st i>ks{(ikki) = 

and we obtain 

^(o) = t y Y 

Itlqlslq i=l dim(y)=t 



ife" 

s — k 



t - k 
s — k 



^kt{ekki) 



t - k' 

s — k 



e-kti 



e k ti(y)e kti (y) 



s— fe. 



rt— fe] rn— 2kl 
Is-kl a l t-k 1 a k (t-k) 



Y\ X \( e kti,e ku ) = \X\h k - 9 t q -q 

i=l Itlqlslq 



Lemma 4.5. With the above notation, 



(17) ^2w(n,s,t;i)P kst (i)P tst (i) = 5 k/ \X\ 2 h k 



□ 



14 CHRISTINE BACHOC, ALBERTO PASSUELLO, FRANK VALLENTIN 

Proof. We compute £ := Y, y ex Ekstfa, y)Ee )S > >t '(y, z). 

h k h t 

S = ^ ^^e fcsi (x)eHi(y)e£ s ' i (y)e tt / j (2;) 

ygX i=l 3=1 
/l fe /if 

= ^ e ksi (x)ei t >j(z) ( y~] e k ti{y)e U j{y) 



tl ) '-IS J 



=1 j=i yex 
^ ^ eksi(x)eu'j{z)\X\(e ktl 

i=l 3=1 

/if; /if 

^ ^ e fesi (x)efl/j(z)|X|(5fc£(5 is '5jj 

i=l 3=1 

.< ; — //.- 

5ki5ts' \X\ 

&kl&ts'\X\ 



) 

n - 2k' 
t - k 



k(t-k) 



n — 


2k 


t - 


k 


n — 


2 k 


t - 


k 



q k{t - k) E kst ,(x,z). 



We obtain, with t = s' , t' = s, x = z & X s , and taking Ea s (y, x) = E( st (x, y) 
into account, 

In - 2k] 



^2 E kst(x,y)E £st (x,y) = 5 M \X\ 
y&Xt 



t - k 



q k(t ~ k) E k 



ss \Xj X) 



The above identity becomes in terms of P ks t 

^2 p kst(s-dim(xr\y))P est (s-dim(xny)) = 5 M \X\ 

yex t 

Now we obtain (|17p by (|16p and Lemma 14.31 



n - 2k 
t - k 



q k{t - k) Pkss(0). 



□ 



We showed that the functions P ks t satisfy the same orthogonality relations 
as the g-Hahn polynomials. So we are done if P ks t is a polynomial of degree 
at most k in the variable [u] = [dim(x n y)]. This property is proved in 
the case s = t in [7, Theorem 5] and extends to s < t with a similar line 
of reasoning. The multiplicative factor between P kst (u) and Q k (n, s,t;u) is 
then given by P ks t(0) and the proof of Theorem 14.21 is completed. □ 

4.2. Symmetry reduction of the semidefinite program (llOj) for pro- 
jective codes. Clearly, (110p is G-invariant: this means that for every fea- 
sible solution F and for every g £ G, also gF is feasible with the same 
objective value. Hence, we can average every feasible solution over G. In 
particular, the optimal value of (I10p is attained by a function F which is 
G-invariant and so we can restrict the optimization variable in (110p to be a 
G-invariant function. 
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A function F(x,y) € R^ x ^ is G-invariant if it depends only on dim(x), 
dim(y), and dim(x n y). So we introduce F, such that F(x,y) = F(s,t,i) 
for x, y E X with dim(x) = s, dim(y) = t, dim(x fl y) = i. Let 

Nsti := \{(x, y) € X x X : dim(x) = s, dim(y) = t, dim(x fl y) = i}\ 

and 

(18) n(d) := {(s,t,i) : < s,t < n, < % < min(s,t), s + t < n + i, 

either s = t = i or s + t — 2i > d}. 

Then, (fTOl) becomes: 

A q (n,d) < max | ^ N sti F(s,t,i) : F e IRN 3 ; p y o, F > 0, 

s.t.i 

n 

^N sss F{s,s,s) = 1, 

F(s,t,i) = if (s,M) 0(d) }, 

where F y means that the corresponding F is positive semidefinite. 
Then, W6 introduce the variables x s ti 

:= N st iF(s,t,i). It is straight- 
forward to rewrite the program in terms of these variables, except for the 
condition F >z 0. From Theorem 14. 1^ this is equivalent to the semidefinite 
conditions F^ y 0, where the matrices F}. are given by the scalar product of 
F and E}.: 

(F k )st= , rn-2fci fe(fl _ fc) rn-2fci fc(t _ fc) E F(x , y)E kst (x , y) 

1 1 K L s— fc J L t— k J <^ u,v,i 

We can substitute the value of Ekst(u,v,i) using Theorem 14. 2\ in particular 
it vanishes when (u,v) ^ {s,t), and, when (u,v) = (s,i) and s < t: 

rt-fc-i 

i Lt J g LsJ q L s— fc J g 

Theorem 4.6. 

^(n,d)<max| ^ x sti : (x s tt)( s ,t,i) e fi(d), x sti > 0, 

^ ^ — 1 ) 

s=0 

i^fc ^ for all k = 0, ... , [n/2\ } , 
where Q(d) is defined in ()18|) and f/ie matrices F k are given in U9\) . 
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Remark 4.7. A projective code C with minimal distance d provides a feasible 
solution of the above program, given by: 

Xsti = j^j \{(x,y) € C : dim(x) = s,dim(y) = t,dim(x n y) = i}. 

In particular, we have 

^Xsti = \C n Q g (n, s)\, 

t,i 

so, we can add the valid inequality 

J2 X °K ^ A q (n,s,2\d/2]) 

t,i 

to the semidefinite program of Theorem \4-6\ in order to tighten it. 

Following the same line of reasoning, we could also add the linear inequal- 
ities 

n r -i 

n 



k 



, k = 0, . . . , n 
q 



^2c(s,k,e)(^2x sti ) < 

s=0 t,i 

where e = [(d — l)/2j , to the semidefinite program of Theorem \4-6\ so that 
the resulting semidefinite program contains all the constraints of the linear 
program of Theorem \2.3l It turns out that this semidefinite program behaves 
numerically badly, and that, when it can be computed, its optimal value 
is equal to the minimum of the optimal values of the initial semidefinite 
program and of the linear program. 

5. Numerical results 

In this section we report the numerical results obtained for the binary 
case q = 2. Table 1 contains upper bounds for A2(n, d) for the subspace dis- 
tance d$ while Table 2 contains upper bounds for A™* 1 (n, d) for the injection 
distance di recently introduced in |17j . 

5.1. Subspace distance. The first column of Table 1 displays the upper 
bound obtained from Etzion-Vardy's linear program, Theorem 12.31 Observ- 
ing that the variables in this program represent integers, its optimal value 
as an integer program gives an upper bound for A q (n, 2e + 1) that may im- 
prove on the optimal value of the linear program in real variables. However, 
we observed a difference with the optimal value of the linear program in real 
variables of at most 1. In Table 1, we display the bound obtained with the 
optimal value of the linear program in real variables, and indicate with a 
superscript * the cases when the integer program gives a better bound (of 
one less). 

The second column contains the upper bound from the semidefinite pro- 
gram of Theorem 14.61 strengthened by the inequalities (see Remark 14.7ft : 

x s ti < A2(n, s, 2\d/2~\) for all s = 0, . . . , n. 

t,i 
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In both programs, A 2 (n, k, 25) 
Theorem E3 
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was replaced by the upper bound from 



parameter 


E-V LP 


SDP 


A 2 (4,3) 


6 


6 


A 2 (5,3) 


20 


20 


A 2 (6,3) 


*124 


124 


^2(7,3) 


832 


776 


^2(7,5) 


36 


35 


A 2 (8,3) 


9365 


9268 


A 2 (8,5) 


361 


360 


A 2 {9,3) 


*114387 


107419 


A 2 (9,5) 


*2531 


2485 


A 2 (10,3) 


*2543747 


2532929 


A 2 (10,5) 


*49451 


49394 


A 2 (10,7) 


*1224 


1223 


A 2 (ll,5) 


693240 


660285 


^2(11,7) 


9120 


8990 


^2(12,7) 


323475 


323374 


A 2 (12,9) 


*4488 


4487 


A 2 (13,7) 


4781932 


4691980 


A 2 (13,9) 


*34591 


34306 


A 2 (U,9) 


2334298 


2334086 


A 2 (14,ll) 


*17160 


17159 


A 2 (15,ll) 


*134687 


134095 


A 2 (16,13) 


*67080 


67079 



Table 1. Bounds for the subspace distance 



5.2. Additional inequalities. Etzion and Vardy [10J found additional valid 
inequalities for their linear program in the special case of n = 5 and d = 3. 
With this, they could improve their bound to the exact value A 2 (5,3) = 18. 
In this section we establish analogous inequalities for other parameters (n, d). 

Theorem 5.1. Let C C 7 7 (F™), of minimal subspace distance d, and let 
D k := \C nQ q (n,k)\. Then, if 

d+2 \d/2] + 2<2n<2d + 2 \d/2] + 2, 

we have: 

• D 2n -d-[d/2]-l < i; 

• if D 2n _ d _ W2] _x = 1 then 

qn _ q2n-d-\d/2\-l 
D W 2 ~\ ^ q \ d/2] _ n -d-l ■ 
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Proof. It is clear that Di < 1 for < i < \d/2] . Moreover, for all x, y € C D 
G(n, \d/2\), x ^ y, dim(xny) = 0. We want to show that D 2n -d-\d/2']-i < !• 
Indeed assume by contradiction x ^ y € C H Q(n,2n — d — \d/2~\ — 1), we 

( An -2d -2 \d/2] - 2 < n + dim(x n y) 
\ d<An-2d-2 \d/2] - 2-2 dim(x n y) 

leading to 

f 2dim(xny) > Qn - Ad - A\d/2~] -A (*) 
\ 2 dim(x n y) < An - 3d - 2 \d/2~] - 2 (**) 

To obtain a contradiction, we must have (*) > (**) which is equivalent to 
the hypothesis 2n > d + 2 \d/2] + 2. 

With a similar reasoning, we prove that, for all x € C n ^(n, [d/2]) and 
w £ C (1 G (n, 2n — d — \d/2~\ — 1), dim(x Pi tf) = n — d — 1. Indeed, 

J 2n — d — l<n + dim(x fl w ) 
| d<2n-d-l-2 dim(x n w) 

so 

J dim(x nw)>n-(i-l 
| dim(a; P\w) < n — d — 1/2 

which yields the result. 

Now we assume D 2n -d- \d/2\-\ = 1- Let w E C n G(n, 2n — d — \d/2~\ — 1). 
Let U denote the union of the subspaces x belonging to C n G(n, \d/2~\). We 
have \U\ = 1 + L> rd/2 ] (g rc(/21 - 1) and |W D «;| = 1 + L> rrf /2l (g™-^ 1 - 1). On 
the other hand, \U\(U Hw))\ < \V£\w\, leading to 

D W 2](q {dm - q n - d ~ l ) <q n - g 2 — rf -r^/2i-i . 

□ 

In several cases, adding these inequalities led to a lower optimal value, 
however we found that only in one case other than (n,d) = (5,3), the final 
result, after rounding down to an integer, is improved. It is the case (n, d) = 
(7, 5), where D% < 17 and, by Theorem l5.lt if D§ = 1 then D% < 16. So we 
can add D 3 + D 5 < 17 and D 2 + D 4 < 17, leading to: A 2 (7, 5) < 34. This 
bound can be obtained with both the linear program of Theorem 12.31 and 
the semidefinite program of Theorem 14.61 

5.3. Injection distance. Recently, a new metric has been considered in 
the framework of projective codes, the injection metric, introduced in |17j . 
The injection distance between two subspaces U, V G V(W™) is defined by 

di(U,V) = max{dim(C/),dim(y)} -dim([/ny). 

When restricted to the Grassmann space, i.e. when U, V have the same 
dimension, the new distance coincides with the subspace distance (up to 
multiplication by 2). In general we have the relation (see [IT] ) 

di(U,V) = \d s {U,V) + \\&m{U) -dim(y)|, 
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where ds denotes the subspace distance. 

It is straightforward to modify the programs in order to produce bounds 
for codes on this new metric space ("P(F"), d,{). Let 

A\ nj (n,d) = max{|C| : C C "P(F"), d l {C) > d}. 

For constant dimension codes, we have A q n:> (n,k,d) = A q (n,k,2d). 

To modify the linear program of Etzion and Vardy for this new distance, 
we need to write down packing-constraints. The cardinality of balls in "P(F™) 
for the injection distance can be found in [14]. Let B tnj (V,e) be the ball 
with center V and radius e. If dim(l/) = i, we have 



\B inj {V,e) 



E 

r=0 

+ EE 



i 




n — i 




r 




r 






q 






r 








— "S 


T 

Q 


(r—a) 


: 


= 1 









i 




n — i 


+ 


i 




n — i 


.) 


r 




r — a 




r — a 




r 






q 








q 







We define c mj (i,k,e) := \B inj (V, e) n Q q (n, k)\ where dim(V) = i. We set 
a := \i — k\. 



c mj (i,k,e) 



r(r—a) ["«] I""- - *] 
r=0 y LrJ q Ir— a\ q 



E 

E e r=0l rir - a) [ r - a ] o [ n 



qL r Iq 



iii>k 
if i < k 



Theorem 5.2 (Linear programming bound for codes in V(¥ q ) with injection 
distance). 

A q nj (n, d) < max | Efc=o xu ■ Xk < A l q nj (n, k, d) V k = 0, . . . , n 

ELo^M^ K],Vfc = 0,..., n } 

For the semidefinite programming bound, we only need to change the 
definition of Vt{d) by 

(20) n inj (d) := {(s,t,i) : < s,t < n, i< mm(s,t), s + t<n + i, 

either s = t = i or max(s, t) — i > d}. 

Theorem 5.3. 

A q nj (n,d) < max | ^ x sti : (x s ti)( SiM)en i»j (d) , a sti > 0, 

(s,(,i)eJ2™J((i) ™ 

/ ^sss — 1 ! 

F fe ^ for allk = 0,..., [n/2\} 

where £l m:, (d) is defined in (120p and i/ie matrices Fk are given in < fi9j) . 

Table 2 displays the numerical computations we obtained from the two 
programs. 
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parameter 


E-V LP 


SDP 


4 nj (7,3) 


37 


37 


4 ni (8,3) 


362 


364 


4 ni (9,3) 


2533 


2536 


4 nj (10,3) 


49586 


49588 


4 ni (10,4) 


1229 


1228 


4 ni (ll,4) 


9124 


9126 


4 nj (12,4) 


323778 


323780 


4 ni (12,5) 


4492 


4492 


4 nj (13,5) 


34596 


34600 


4^(14,6) 


17167 


17164 


# 2 nj (15,6) 


134694 


134698 


4 nj (16,7) 


67087 


67084 



Table 2. Bounds for the injection distance 



Remark 5.4. We observe that the bound obtained for A 2 (n, 2e + l) is most 
of the time slightly larger than the one obtained for A2(n,4e + 1). In [13], the 
authors noticed that their constructions led to codes that are slightly better for 
the injection distance that for the subspace distance. So both experimental 
observations indicate that A2(n, 4e + 1) is larger than A™ 3 (n, 2e + 1). 

The computational part of this research would not have been possible 
without the use of free software: We computed the values of the linear pro- 
grams with Avis' Irs 4.2 available from http://cgm.cs.mcgill.ca/~avis/C/lrs.html. 
The values of the semidefinite programs we computed with SDPA or SDPT3, 
available from the NEOS website http://www.neos-server.org/neos/. 
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